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本 振 相 位 噪声 及 其 对 CMSK 解 调 性 能 的 影响 
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摘要 :GMSK 调制 方式 在 空间 通信 中 的 应 用 越 来 越 受 重视 ,本 振 相 位 噪声 对 GMSK 解 调 具有 不 可 忽略 的 
影响 。 在 分 析 相 位 噪声 产生 原理 的 基础 上 ,推导 了 其 时 域 和 频 域 表达 式 , 并 对 其 进行 了 仿真 验证 。 通 过 
在 GMSK 的 Simulink 调制 解 调 系 统 中 加 入 相位 噪声 ,仿真 分 析 了 不 同 程度 的 相位 骂 声 对 解 调 性 能 的 影 
响 。 研 究 方 法 和 结论 可 以 推广 到 存在 相位 噪声 的 任意 调制 类 型 的 解 调 系 统 ,对 解 调 系 统 的 优化 具有 重 
要 参考 价值 。 
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Analysis of Local Oscillator Phase Noise and Its Effect 
On the Demodulation Performance of GMSK 
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Abstract: GMSK modulation plays an important role in space telecommunications. Local osciliator phase noise 
has significant influence on GMSK demodulation system. Some characteristics of phase noise and its math ex- 
pressions on both time and frequency domains are analyzed. Furthermore, the simulated phase noise sequence 
is added to the GMSK modulation system to show how the different levels of phase noise affect the demodula- 


tion performance. The conclusions have deep reference value to optimize the GMSK demodulation system. This 


method to analyze phase noise and its effect can be extended to any modulation system. 
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在 通信 系统 中 ,一 个 理想 的 本 振 输出 信和 号 频谱 是 
一 根 单 谱 线 ;而 实际 的 本 振 不 可 避免 地 存在 随机 相位 
持 动 ( 即 相 位 噪声 ) ,通信 系统 的 误 码 率 将 恶化 。 提 高 
本 振 的 频谱 纯度 可 以 改善 通信 系统 的 误 码 性 能 。 而 当 
误 码 率 达 到 一 定数 量 级 时 ,过 分 追求 高 性 能 的 本 振 并 
不 能 显著 提高 系统 误 码 性 能 ,反而 会 使 本 振 的 造价 变 
得 非常 昂贵 甚至 不 可 实现 "。 因 此 ,分 析 通 信 系统 中 
本 振 相 位 噪声 对 误 码 性 能 的 影响 ,有 助 于 给 工程 设计 
提供 有 效 参考 ,以 便 兼顾 系统 的 可 实现 性 和 高 性 能 。 

GMSK 调制 方式 具有 恒 包 络 的 特性 ,适用 于 饱和 
功放 , 且 带 宽 效 率 高 , 因此 GMSK 调制 方式 已 经 被 
CCSDS 建议 为 深 空 通信 高速 数 传 的 首选 调制 方式 。 
ESA 和 JAXA 合作 研发 的 计划 在 2014 年 发 射 的 水 星 探 


收 稿 日 期 :2014 -01 -05 
方向 为 空 闻 通信 系统 。 


测 需 贝 皮 -哥伦布 (bepicolombo) 已 经 采用 GMSK 调制 
体制 ,这 是 在 深 空 应 答 机 上 首次 使 用 CMSK 调制 方式 ， 
同时 国内 对 空间 通信 中 的 GMSK 的 研究 也 越 来 越 多 ” 。 

笔者 将 分 析 相 位 噪声 的 原理 ,在 此 基础 上 推导 其 
时 域 和 频 域 表 达 式 ,并 进行 仿真 分 析 。 通 过 搭建 GM- 
SK 调制 解 调 系统 的 Simulink 模型 ,分析 不 同 程度 的 相 
位 噪声 对 解 调 系统 的 影响 。 
1 相位 噪声 原理 分 析 

一 个 理想 的 本 振 输 出 信号 可 以 表示 为 

V(t) = Vsin(2mh1t) (1) 

其 中 所 均 为 常数 。 这 是 一 个 纯净 的 本 振 信 号 ,从 
频 域 看 本 振 信 号 为 一 根 纯净 谱 线 ,从 时 域 看 正 牙 波 的 
周期 为 一 恒定 值 。 但 实际 中 任何 信号 不 可 避免 地 被 噪 
声调 制 。 

实际 本 振 可 表示 为 


Vt) =[ WV +e(z) Jsin[L 2mfos + p(t)] +V (7) (2) 
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9(t) 为 各 种 原因 造成 的 附加 瞬时 相位 起 伏 ,通常 情况 
下 (zt) << WW,Vi(i) << WW, 即 在 接近 所 的 偏 移 频率 
上 , 式 (2) 可 以 简化 为 


V(t) = Wsin[2m6t+e(D] (3) 
本 振 输 出 的 瞬时 频率 可 以 记 为 
ot) = 及 + 于 电 全 (4) 
相对 瞬时 频率 抖动 可 定义 为 
1 = (5) 
相对 瞬时 相位 拌 动 可 定义 为 


(1) =2 (6) 


由 晶振 自身 引起 的 随机 相位 抖动 就 是 相位 噪声 。 
y(t) 包含 了 频率 稳定 度 的 全 部 信息 ,参考 文献 
[3] ~ 文献 [5] 认 为 ,y(i) 的 功率 谱 密 度 5,(f) 可 以 由 
香 律 谱 噪声 模型 表示 。 
传统 的 宪 律 模型 如 下 所 示 : 
SO = Dhf’ (7) 
式 中 ,hh 为 噪声 指数 是 的 强度 系数 。8 取 值 从 ( -2， 
-1,0,1,2) 变 化 Be] 。 虽 然 在 传统 的 宕 律 模型 中 ,只 考 
虑 了 取 整 数值 的 情况 ,但 是 在 广泛 的 实验 验证 中 , 
为 非 整 数值 的 噪声 源 也 不 容 忽 视 ']。 
由 于 x(?) 是 y(t) 的 积分 ,所 以 x(i) 的 功率 谱 密 度 
与 5,(f) 的 关系 为 
AR 
3,.(f) 二 (2m 站 : 
为 了 简化 公式 ,重新 定义 变量 a, 可 以 得 到 


I h, b-2 oc a 
$5. -Ta “1 (9) 


本 振 相 位 曝 声 具有 1/ 挛 形式 ,通常 考虑 其 范围 为 0 < 
Q <2。 

由 于 这 些 寡 律 型 噪声 基本 不 满足 平稳 遍历 条 件 ， 
它们 的 标准 方差 是 不 收敛 的 ,因而 其 频率 稳定 度 的 时 
域 表征 需要 通过 Allan 方差 9 修正 Allan 方差 ""、 
Hadamard 方差 等 来 表征 '*1。 


2 相位 噪声 序列 生成 


2.1 分 形 布朗 运动 模型 

基于 多 年 以 来 相关 学 者 的 大 量 研究 ,常用 的 产生 
1/f “噪声 的 连续 建 模 方 法 被 分 为 三 大 类 "] :ARMA 模 
型 ,分 形 积分 ,小 波 分 析 。 本 文 的 离散 模型 是 基于 分 形 
积分 的 方法 。 

分 形 之 父 Mandelbrot 观察 到 1 "噪声 和 布朗 运 
动 的 功率 谱 密 度 存在 一 定 的 相似 性 。 布 朗 运动 的 功率 
谱 密 度 是 正比 于 1/f? 的 , Mandelbrot 将 布朗 运动 的 概 
念 进 行 推广 ,将 1/f “噪声 看 做 一 种 分 形 布朗 运动 


(8) 
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(FBM ,fractional Brownian motion) "| 。 

布朗 运动 B(i) 是 零 均 值 的 随机 过 程 , 且 满 足以 下 
3 个 性 质 : 

中 8(0) =0; 

@) B(1) -B(s) 是 平稳 的 , 且 均 值 为 0, 方差 正比 
于 (1 -$s) ,Ss<t; 

Bs) - B(H),B(is) -Bh),, Bi) 一 
B(t,_1) 是 独立 的 ,ti 和 和 …t。 

分 形 布 朗 运动 B.(it) 可 由 此 定义 为 满足 以 下 3 个 
性 质 的 零 均值 随机 过 程 : 

由 Br(0) =0; 

@ B,(t) -Bj(s) 是 平稳 的 , 且 均 值 为 0, 方 差 正比 
于 (1-s)"! ,set; 

® Bi,) -Bi ) ,Br(ts) -Bt,) (人 > 
Br 1) 是 独立 的 ,bi 志和; 

从 大 量 的 实验 和 观察 证 明 出 1/ 严 噪声 是 自然 界 的 
基本 性 质 ,是 普遍 存在 的 , Mandelbrot 将 1/ 户 噪声 归 类 
为 一 种 随机 分 形 ,因此 ,还 加 了 第 4 条 性 质 ; 

@ Bj(t) 是 自 相似 的 (或 者 说 标 度 不 变 ) ,其 标 度 
不 变 是 针对 功率 谱 密 度 来 说 的 ,假如 频率 的 标 度 变 了 ， 
原先 的 幅度 标 度 可 以 由 直接 乘 以 一 个 常数 来 得 到 。 
2.2 ”相位 噪声 模型 推导 
2.2.1 和 白 嗓 声 与 布衣 运动 

由 文献 [5] 可 知 ,布朗 运动 B(it) 是 白 噪声 w(t) 的 
积分 ,其 表达 式 为 

Blt) = fea (10) 
对 式 (10) 进 行 Laplace 变换 得 到 
H(s) = 二 (11) 
其 对 应 的 离散 差分 方程 为 
B, = B，， +w, (12) 
对 式 (12) 做 z 变换 ,得 到 传递 函数 


H(z) (13) 


2.2.2 相位 骂 声 与 白 品 声 
布朗 运动 对 应 到 式 (9) 中 即 是 a =2 的 情况 。 正 
如 白 噪声 是 布朗 运动 的 微分 ,相位 噪声 x( 电 可 由 以 下 
分 数 阶 微分 来 定义 中 
xl) =w(1) (14) 


其 中 wa 是 实数 ,范围 为 0<aw<2。 由 黎 曼 - 刘 维 尔 分 
形 积分 公式 可 以 得 到 式 (14) 的 解 为 
x(£) = pe 7) (7r) dr (15) 


这 等 效 于 一 个 由 白 噪 声 激励 的 线性 系统 ,其 冲 激 
响应 为 
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h(t) -ay (16) 
经 过 Laplace 变换 和 归 一 化 ,得 到 
已 (= (17) 
由 冲 激 不 变法 ,可 得 到 离散 模型 传递 函数 为 
Bs) = 0 re>1) (18) 


离散 的 相位 噪声 序列 x; ,可 以 由 高 斯 白 噪声 序列 
w; 通过 传递 函数 为 瓦 (z) 的 线性 系统 来 表示 ,其 z 域 
及 时 域 表 达 式 为 


X(z) =H(z) W(z) (19) 
NX 二 > (20) 

由 文献 [5] 得 到 式 (19) 对 应 的 功率 谱 密度 为 
Se =AtQH(z) H(z ) (21) 


其 中 ,Q。 是 输入 的 高 斯 白 噪声 序列 w 的 方差 , 乘 以 At 
是 为 了 使 连续 系统 和 离散 系统 的 标 度 一 致 。 

将 式 (18) 和 z=e”w 代 入 式 (21) ,得 到 功率 谱 密 
度 为 


OA 
$5. = (2sin( nfAL) ) 20 
对 于 远 小 于 奈 奎 斯 特 频 率 的 情况 ,可 近似 为 
OA 
5D sp (23) 
把 式 (18) 按 级 数 展开 ,有 
H(z) =1 | + + + (24) 
式 (24) 中 ,各 项 系数 可 以 由 如 下 递归 算法 得 到 : 
加 = 有 =( 全 + (25) 


由 此 ,采用 式 (25) 可 计算 出 冲击 响应 ,代入 式 
(20) 可 以 直接 产生 相位 噪声 序列 。 


3 相位 噪声 仿真 分 析 


1/F 形式 的 相位 噪声 , 取 a = 1, 在 Simulink 中 搭 
建 相应 模型 ,由 高 斯 白 噪 声 通 过 FIR 滤波 器 ,生成 1 
噪声 , 其 中 滤波 器 的 系数 由 式 (25 ) 得 到 ,并 且 除 以 一 
个 常数 ,使 得 由 滤波 器 引入 的 功率 增益 为 1, 以便 控制 
高 斯 噪声 的 0 就 能 直接 控制 相位 噪声 的 Allan 方差 。 
生成 的 噪声 序列 如 图 1 所 示 ,长 度 为 10000。 

求 出 该 序列 的 功率 谱 密度 ,与 式 (23 ) 的 理论 值 进 
行 比较 ,如 图 2 所 示 。 

可 以 看 到 功率 谱 密度 仿真 值 和 理论 值 接近 ,斜率 
都 为 1 ,符合 If 特性 。 

然而 在 实际 的 系统 中 ,相位 噪声 的 频谱 并 非 在 全 
频率 范围 内 都 是 一 个 固定 的 值 ,其 频率 特性 通常 是 近 
似 分 段 线性 的 ,以 国家 军用 标准 为 例 , 其 参数 如 表 1 所 
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示 , 其 频谱 被 分 为 4 段 。 


0 5000 10000 


图 1 wa=1 时 的 噪声 序列 


一 一 PDS simulation 
+ PDS theory 


图 2 功率 谱 密度 5,(f) 


表 1 国家 军用 标准 相 曲 参数 


偏 移 载 频频 率 /kHz LL( 有 )/dBe:Hz 频谱 分 段 斜率 
0.1 -55 1.5 
1 -70 0.7 
10 -77 0.5 
100 -82 1.8 
1000 -100 
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图 3 用 Simulink 产生 分 段 线性 的 相位 噪声 


图 3 中 的 频谱 仪 可 以 显示 用 该 方法 产生 的 相位 品 
声 的 频谱 图 。 如 图 4 所 示 。 
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图 4 相位 噪声 谱 
将 图 4 和 表 1 进行 对 比 ,可 看 出 仿真 结果 和 理论 
值 接近 ,验证 了 该 方法 的 可 行 性 。 


4 相位 噪声 对 GMSK 解 调 性 能 的 影响 


对 于 空间 通信 系统 ,要 求解 调 器 具有 最 小 的 解 调 
损失 ,并 且 能 够 提取 载波 ,以 此 来 测量 积分 多 普 勒 从 而 
用 于 导航 "" 1。 因此 ,空间 通信 中 GMSK 调制 系统 需要 
采用 相干 解 调 方案 。 

在 相干 解 调 的 Simulink 仿真 中 ,高 斯 信道 模块 品 
声 模式 ELNO 取 0 ~ 11 , 步 进 值 为 1, 误 码 率 统计 模块 
可 输出 在 各 个 E65/N0 取 值 所 对 应 的 误 码 率 结果 变量 。 
图 3 生成 的 相位 噪声 加 在 发 端 NCO 上 。 

GMSK 系统 主要 参数 设置 如 表 2 所 示 。 

表 2 GMSK 系统 主要 参数 


码 速率 /Mbit*s-! 中 心 频率 /MHz 采样 率 /MHz 
1 5 100 


100 200 
Frame:310 


© 


取 6 组 不 同 的 数据 进行 实验 ,如 表 3 所 示 。 分 别 
仿真 无 相位 噪声 、 相 品 小 于 国 军 标 (实验 一 .实验 二 ) 、 
相 品 恰好 满足 国 军 标 ( 实验 三 )、 相 品 大 于 国 军 标 ( 实 
验 四 、 实 验 五 ) 的 情况 下 对 GMSK 误 码 性 能 的 影响 。 
其 中 均 方 根 拌 动 与 功率 谱 密 度 L(/) 的 关系 可 参考 文 
献 [12]。 

表 3 Simulink 仿真 相 品 参数 设置 


L(f) /dBe: Hz ~! 均 方 根 拌 

偏 移 载 频频 率 /Hz 100 1lk 10k 100k 1M 动 /(") 

元 相 噪 NULL 0 
实验 一 -61 -76 -83 -88 -106 1.94 

有 实验 二 -58 -73 -80 -85 -103 2.75 
相 实验 三 ( 国 军 标 ) -55 -70 -77 -82 -100 3.88 
品 实验 四 -52 -67 -74 -79 -97 5.48 
实验 五 -49 -64 -7?71 -76 -94 7.73 

各 实验 误 码 率 统计 结果 如 图 5 所 示 。 


由 图 5 可 以 看 出 , 当 码 速率 为 1 Mbit/s 时 ,相位 噪 
声 对 解 调 性 能 有 较 大 影响 ,可 得 出 如 下 明显 结论 : 


。0S ， 


Error Rate 


Eb/NO/dB 


图 5 误 码 率 统计 


Q@ 即便 相 品 功率 谱 密度 满足 国 军 标 ,时 域 的 均 方 
根 拉动 也 达到 了 3. 88。, 误 码 率 在 10 处 的 解 调 损耗 
已 经 大 于 1 dB。 

@) 实验 一 .实验 二 的 解 调 损耗 都 小 于 1 dB ,结合 
表 3 可 以 看 出 ,实验 一 .实验 二 分 别 对 应 于 相 品 功率 谱 
密度 比 国 军 标 小 6 dB 和 3 dB。 因 此 ,在 工程 应 用 上 ，， 
应 尽量 将 偏离 载 频 1 MHz 以 内 的 相位 噪声 控制 在 至 
少 比 国 军 标 低 3 dB 的 范围 以 下 ,以 保证 解 调 损耗 小 于 
1 dB。 

@ 从 实验 四 、 实 验 五 可 以 看 出 , 随 着 相位 噪声 的 
继续 增 大 ,BER 曲线 严重 恶化 ,提高 到 /NO 值 对 于 系 
统 误 码 率 的 改善 也 很 小 , 误 码 率 出 现 了 不 归 零 的 趋势 。 


5 结束语 


在 分 析 了 相位 噪声 的 基础 上 , 对 相位 噪声 序列 进 

行 仿真 ,并 在 Simulink 平台 仿真 不 同 程度 的 相位 噪声 

对 GMSK 解 调 性 能 的 影响 ,可 以 得 出 如 下 结论 :如 果 

要 使 相位 噪声 对 解 调 系统 的 影响 很 小 , 需 保 持 相 噪 功 

率 谱 密度 比 国 军 标 低 3 dB 的 范围 以 下 。 

用 本 文 的 分 析 和 仿真 方法 ,对 相位 噪声 对 解 调 损 

耗 的 影响 做 前 期 的 定量 仿真 分 析 , 可 以 为 系统 设计 环 

节 中 相位 噪声 的 指标 分 配 提 供 一 定 的 参考 依据 。 该 方 

法 不 仅 适 用 于 GMSK 调制 解 调 系统 ,而且 可 以 推广 到 

存在 相位 噪声 的 任意 调制 类 型 的 解 调 ,具有 普遍 意义 。 
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是 否 初次 使 用 ? 四 
| 使 用 上 次 的 设置 
达到 触发 状态 ? 


图 3 人 机 交互 模块 程序 流程 图 


3 系统 测试 分 析 及 改进 


本 系统 基于 上 文 原理 完成 实物 并 进行 多 次 测试 。 
测试 场景 选择 在 放置 了 一 台 柜 式 格力 空调 的 普通 房间 
内 ,在 GSM 信号 接收 正常 及 红外 发 射 头 正 对 空调 红外 
接收 区 的 情况 下 ,分 别 发 送 温度 从 16 ~30 的 制冷 及 
制 热 模 式 开启 空调 的 指令 ,指令 完成 率 为 100% 且 设 
备 距 空调 可 达 数 米 ; 但 在 红外 发 射 头 与 空调 接收 区 呈 
现 一 定 的 角度 或 其 之 间 有 小 部 分 的 障碍 物 时 ,指令 完 
成 率 有 所 下 降 。 

对 多 次 测试 后 的 结果 分 析 , 仅 有 单个 红外 发 射 头 
及 较 差 的 网 络 环境 使 指令 完成 率 不 能 达到 用 户 的 期 户 
要 求 , 故 将 红外 发 射 模块 及 TC35 的 天 线 部 分 进行 了 
改进 。 在 红外 发 射 模块 方面 ,根据 每 个 发 射 头 所 能 覆 
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盖 到 的 有 效 角度 范围 进行 建 模 ,在 设备 的 重要 面 上 加 
装 了 5 个 红外 发 射 头 ,同时 小 幅度 改进 了 驱动 电路 ,使 
红外 信号 可 覆盖 整个 房间 且 辐 射 功 率 小 幅 上 升 ;在 网 
络 信号 方面 ,改进 了 设备 样式 ,将 天 线 加 强 后 伸 出 设 
备 ,使 网 络 信号 大 幅 优化 。 系 统 改进 后 ,设备 的 放置 位 
置 更 加 灵活 ,设备 与 空调 的 有 效 距离 有 所 提升 ,出 现 网 
络 错误 的 概率 减 小 。 最 后 在 放置 挂 式 或 柜 式 空调 的 房 
间 内 均 进行 测试 ,指令 完成 率 接近 100% 。 


4 结束 语 


智能 远程 家 居 控 制 系统 的 设计 使 在 任何 有 GSM 

网 络 的 地 方 都 能 对 家 居 设 备 进 行 监控 ,加 快 了 家 庭 信 

息 化 、 自 动 化 的 进程 。 利 用 上 文 的 思路 可 以 继续 添加 

厨房 火警 煤气 浓度 检测 、 防 盗 等 模块 ,同时 本 系统 也 

可 以 向 工厂 \ 学 校 等 地 方 扩 展 , 有 很 高 的 可 扩展 性 , 同 

时 突破 了 传统 的 有 线 控制 ,具有 成 本 低 、 智 能 化 高 的 特 

点 ,有 很 高 的 市 场 利用 价值 。 
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